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We compare a set of equations of state derived within microscopic many-body approaches, and study their 
predictions as far as phenomenological data on nuclei from heavy ion collisions, and astrophysical observations 
on neutron stars are concerned. We find that all the data, taken together, put strong constraints not easy to 
be fulfilled accurately. However, no major discrepancies are found among the selected equations of state and 
with respect to the data. The results provide an estimate of the uncertainty on the theoretical prediction at a 
microscopic level of the nuclear equation of state. 
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I. INTRODUCTION 



A convergent effort of experimental and theoretical nuclear 
physics has been developing along several years to determine 
or to constrain the Equation of State (EoS) of nuclear matter. 
Two main areas of research have provided relevant hints in this 
direction. The experimental data on heavy ion collisions have 
been systematically analyzed on the basis of detailed simula- 
tions, in order to constrain not only the density dependence of 
the EoS, but also its isospin dependence. Excellent reviews 
on this subject can be found in the literature [0, 01. On the 
other hand, the analysis of astrophysical observational data, 
noticeably on compact objects, are of great relevance for the 
study of the nuclear EoS. The results in this field of research 
can be considered complementary to the ones that can be ob- 
tained within the heavy ion research activity, because nuclear 
matter is involved in different physical conditions. In Neu- 
tron Stars (NS) nuclear matter is present in beta equilibrium 
from very low density to several times saturation density, and 
it is therefore extremely asymmetric, much more than nuclei 
in laboratory. Despite the different physical situations, an ac- 
curate microscopic theory of nuclear matter is expected to be 
able to explain correctly the data obtained in both physical 
realms. For a pedagogical introduction to the field see ref. 
|0]. Another more phenomenological approach is based on 
the Energy Density Fuctional (EDF) method, in particular the 
Skyrme forces scheme. In this approach a phenomenological 
force or EDF includes a certain set of parameters that are fixed 
by fitting the binding energy of nuclei throughout the mass ta- 
ble. Some of these forces are adjusted also to microscopic 
nuclear matter EoS. Recently, an ample set of Skyrme forces, 
that have been presented in the literature, has been analyzed 
and confronted with the available constraints on the nuclear 
EoS obtained from heavy ion reactions and astrophysical ob- 
jects (j]. The few Skyrme forces that passed these tests have 
been then used to predict the binding energy of a wide set of 
nuclei [5], and none seems to perform satisfactorily well in 
this case. These results show clearly that it is not at all trivial 
to satisfy the constraints coming from astrophysical and heavy 
ion data, as it was already found in ref.fl. Furthermore, if the 
nuclear mass data are included, it seems very difficult to re- 
produce all the data sets. It looks that the phenomenological 
forces are not flexible enough to this purpose. In this paper we 



follow a complementary approach. We focus on the micro- 
scopic theory of nuclear matter and the corresponding EoS's 
that have been developed so far. Adopting the same line as in 
ref. 01, we test the predictions of different microscopic many- 
body approaches with respect to the constraints coming from 
experiments and phenomenology. The aim of this analysis is 
to establish to what extent the data can be reproduced within a 
microscopic approach and the accuracy of the present theory 
of nuclear matter. 

This paper is organized as follows. In Sect. II a review of the 
currently used microscopic many-body methods is presented, 
whereas results are discussed in Sect. III. Finally, in Sect. IV 
we draw our conclusions. 



H. EOS OF NUCLEAR MATTER IN MICROSCOPIC 
APPROACHES 



Empirical properties of infinite nuclear matter can be cal- 
culated using many different theoretical approaches. In this 
paper we concentrate on the microscopic ones, the only in- 
put being a realistic free nucleon-nucleon (NN) interaction 
with parameters fitted to NN scattering phase shifts in differ- 
ent partial wave channels, and to properties of the deuteron. 
In the following we discuss the non-relativistic Brueckner- 
Hartree-Fock (BHF) method [7] and its relativistic counter- 
part, the Dirac-Bruckner- Hartree-Fock (DBHF) approxima- 
tion QSJ] and the Variational method H. A relatively more 
recent approach is based on the Chiral effective Field Theory, 
which has been extensively applied to Nuclear Matter. For a 
review see ref. illoL fllTl - It establishes a link with the under- 
lying QCD structure of strong interaction, and it has reached 
a high degree of sophistication ifl^ - THill . We will not discuss 
this approach, because the corresponding EoS is based on a 
low momentum expansion, and therefore its behaviour at high 
density is uncertain. Hence, a comparison with the other mi- 
croscopic EoS's would be incomplete. However, the EoS from 
the Chiral approach turns out to be close to the variational EoS 
lHoll . and therefore the parameters near saturation are expected 
to be quite similar. 
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A. The Brueckner-Bethe-Goldstone approach 

The Brueckner-Bethe-Goldstone (BBG) theory is based on 
a linked cluster expansion of the energy per nucleon of nuclear 
matter (see Ref. [7], chapter 1 and references therein). The 
basic ingredient in this many-body approach is the Brueck- 
ner reaction matrix G, which is the solution of the Bethe- 
Goldstone equation 

/-r„ ^i , V \kah)Q{k a k b \ 
G[p;co]=v + j: V G[p;co], (1) 

where v is the bare NN interaction, p is the nucleon number 
density, and co the starting energy. The single-particle energy 
e(k) (assuming h=\ here and throughout the paper), 

k 2 

e(k)=e(k;p) = — + U(k;p), (2) 

and the Pauli operator Q determine the propagation of inter- 
mediate baryon pairs. The Brueckner-Hartree-Fock (BHF) 
approximation for the single-particle potential U{k\p) using 
the continuous choice is 

U(k;p)=Re £ {kk!\G\p;e{k)+e{k!)]\kk!) a , (3) 

k'<k F 

where the subscript "a" indicates antisymmetrization of the 
matrix element. Due to the occurrence of U{k) in Eq. ©, 
they constitute a coupled system that has to be solved in a 
self-consistent manner for several momenta of the particles 
involved, at the considered densities. In the BHF approxima- 
tion the energy per nucleon is 

7 = \^ + h £ (kk'lGlp-X^+eik'mkX. (4) 
A 5 2m 2p k ^ kp 

In this scheme, the only input quantity we need is the bare NN 
interaction v in the Bethe-Goldstone equation (1). The nuclear 
EoS can be calculated with good accuracy in the Brueckner 
two hole-line approximation with the continuous choice for 
the single-particle potential, since the results in this scheme 
are quite close to the calculations which include also the three 
hole-line contribution II 1 711 . 

However, it is commonly known that non-relativistic calcu- 
lations, based on purely two-body interactions, fail to repro- 
duce the correct saturation point of symmetric nuclear matter. 
One of the well known results of several studies, that lasted 
for about half a century, is the need of introducing three-body 
forces (TBFs). In our approach the TBF is reduced to a den- 
sity dependent two-body force by averaging over the position 
of the third particle, assuming that the probability of having 
two particles at a given distance is reduced according to the 
two-body correlation function [18]. 

In this work we will illustrate results for two different ap- 
proaches to the TBF's, i.e. a phenomenological and a micro- 
scopic one. The phenomenological approach is based on the 
so-called Urbana model, which consists of an attractive term 



due to two-pion exchange with excitation of an intermediate 
A resonance, and a repulsive phenomenological central term 
lfl9ll . We introduced the same Urbana three-nucleon model 
within the BHF approach. Those TBF's produce a shift of 
about +1 MeV in energy and —0.01 fm ~ 3 in density. This ad- 
justment is obtained by tuning the two parameters contained 
in the TBF's, and was performed to get an optimal saturation 
point (the minimum) (for details see Ref. IU8IP . 

The connection between two-body and three-body forces 
within the meson-nucleon theory of nuclear interaction is 
extensively discussed and developed in references ll20l l2~fll . 
At present the theoretical status of microscopically derived 
TBF's is still quite rudimentary, however a tentative approach 
has been proposed using the same meson-exchange parame- 
ters as the underlying NN potential. Results have been ob- 
tained with the Argonne vig, the Bonn B, and the Nijmegen 
93 potentials HI. 

In the past years, the BHF approach has been extended in 
order to include the hyperon degrees of freedom [22], which 
play an important role in the study of neutron star matter. 
However, in this paper we are mainly interested on the prop- 
erties of the nucleonic EoS, therefore this issue will not be 
discussed further. 



B. The Relativistic approach 

The relativistic framework is the one on which the nuclear 
EoS should be ultimately based. The best relativistic treat- 
ment developed so far is the Dirac-Brueckner (DBHF) ap- 
proach 1 8]. The DBHF method can be developed in analogy 
with the non-relativistic case, i.e. the nucleon inside the nu- 
clear medium is viewed as a dressed particle in consequence 
of its two-body interaction with the surrounding nucleons. 
The two-body correlations are described by introducing the 
in-medium relativistic G-matrix. The DBHF scheme can be 
formulated as a self-consistent problem between the single 
particle self-energy E and the G-matrix. Schematically, the 
equations can be written 

G = V + iJ VQggG (5) 

E = -iJ F (Tr[gG]~gG) (6) 

where Q is the Pauli operator which projects the intermediate 
two particle momenta outside the Fermi sphere, as in the BBG 
G-matrix equation ([TJ, and g is the single particle Green's 
function, which fulfills the Dyson equation 

g=go+goZg (7) 

where go is the (relativistic) single particle Green's function 
for a free gas of nucleons, and E is the nucleon self-energy 
which expresses the influence of the surrounding nucleons. 
The self-energy can be expanded in the covariant form 

Z(k, kf ) = E,. (Jt, k F ) - 70E0 (k,k F ) + y kE v (8) 

where 7^ are the Dirac gamma matrices, and the coefficients 
of the expansion are scalar functions, which in general depend 
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on the modulus |k| of the three-momentum and on the energy 
Icq. The free single particle eigenstates, which determine the 
spectral representation of the free Green's function, are solu- 
tions of the Dirac equation 



[ Yn^ -M]u{k) = 



(9) 



where u is the Dirac spinor at four-momentum k. For the full 
single particle Green's function g the corresponding eigen- 
states satisfy 



-M + L]u{k)* = 



(10) 



Inserting the above general expression for E, eq.®, after a 
little manipulation one gets 



with 



[ 7^* -M* ]u{kf = 



1+E V 



; k' = k' ; M* 



M + E, 
1+E V 



(11) 



(12) 



This is the Dirac equation for a single particle in the medium, 
and the corresponding solution is the spinor 



u*(k,s) = 



-M* 



2M* 



1 

a ■ k 



6 ; E£ = Vk 2 + M* 



(13) 

In line with the Brueckner scheme, within the BBG expan- 
sion, in the self-energy of equation (O only the contribution 
of the single particle Green's function pole is considered. Fur- 
thermore, negative energy states are neglected and one gets 
the usual self-consistent condition between self-energy and 
scattering G-matrix. 

In any case, the medium effect on the spinor of equa- 
tion (fT3l is to replace the vacuum value of the nucleon mass 
and three-momentum with the in-medium values of equation 
(fT2l . This means that the in-medium Dirac spinor is "ro- 
tated" with respect to the corresponding one in vacuum, and a 
positive (particle) energy state in the medium has some non- 
zero component on the negative (anti-particle) energy state 
in vacuum. In terms of vacuum single nucleon states, the 
nuclear medium produces automatically anti-nucleon states 
which contribute to the self-energy and to the total energy of 
the system. It has been shown in ref. M23I1 that this relativistic 
effect is equivalent to the introduction of well defined TBF's 
at the non-relativistic level. These TBF's turn out to be repul- 
sive, and consequently produce a saturating effect. Actually, 
including in BHF only these particular TBF's, one gets results 
close to DBHF calculations, see ref.|24|. Generally speaking, 
the DBHF gives in general a better saturation point than BHF, 
and the corresponding EoS turns out to be stiffer above satu- 
ration than the one calculated from the BHF + TBF method. 

In the relativistic context the only NN potentials which have 
been developed are the ones of one boson exchange (OBE) 
type. In the calculations shown here the Bonn A potential is 
used HI. 



C. The Variational method 

In the variational method H one assumes that the ground 
state wave function SP can be written in the form 

^i,r 2 , ) =U i< jf(r i j)0(r 1 ,r 2 , ), (14) 

where <t> is the unperturbed ground state wave function, prop- 
erly antisymmetrized, and the product runs over all possible 
distinct pairs of particles. The correlation factor is here deter- 
mined by the variational principle, i.e. by imposing that the 
mean value of the Hamiltonian gets a minimum 
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Sf 



= 0. 



(15) 



In principle this is a functional equation for the correlation 
function /, which however can be written explicitly in a 
closed form only if additional suitable approximations are in- 
troduced. The function f(rij) is assumed to converge to 1 at 
large distance and to go rapidly to zero as ru — > 0, to take into 
account the repulsive hard core of the NN interaction. Further- 
more, at distance just above the core radius a possible increase 
of the correlation function beyond the value 1 is possible. 

For nuclear matter it is necessary to introduce a channel de- 
pendent correlation factor, which is equivalent to assume that 
/ is actually a two-body operator Fjj. One then assumes that 
F can be expanded in the same spin-isospin, spin-orbit and 
tensor operators appearing in the NN interaction. Momentum 
dependent operators, like spin-orbit, are usually treated sepa- 
rately. The product in equation (fl4l must be then symmetrized 
since the different terms do not commute anymore. 

If the two-body NN interaction is local and central, its mean 
value is directly related to the pair distribution function g(r) 



<V>=\p I cF rv(r)g[T) . 



where 



g0i-r 2 ) 



jU i>2 d 3 ri \^( ni r 2 ....)\ 2 
fn#Pr t \W(n,Ti~..)\ 2 



(16) 



(17) 



The main job in the variational method is to relate the pair 
distribution function to the correlation factors F. Again, in nu- 
clear matter also the pair distribution function must be consid- 
ered channel dependent and the relation with the correlation 
factor becomes more complex. In general this relation can- 
not be worked out exactly, and one has to rely on some suit- 
able expansion. Furthermore, three-body or higher correlation 
functions must in general be introduced, which will depend 
on three or more particle coordinates and describe higher or- 
der correlations in the medium. Many excellent review papers 
exist in the literature on the variational method and its exten- 
sive use for the determination of nuclear matter EoS [9, 26]. 
The best known and most used variational nuclear matter EoS 
is the one calculated by Akmal, Pandharipande and Ravenhall 
B27I1 . In their paper the authors showed calculations using the 
Argonne Vi8 NN interaction [28], with boost corrections to 
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the two-nucleon interaction, which give the leading relativis- 
tic effect of order (v/c) 2 , as well as three-nucleon interactions 
modeled with the Urbana force explained above. This EoS 
will be tested in the present paper. 



III. RESULTS AND DISCUSSION 

A. Phenomenology of nuclei 

One of the main goals of the microscopic many -body meth- 
ods is the correct reproduction of the nuclear matter satu- 
ration point. The non-relativistic approaches typically lead 
to an over prediction of the saturation density po = 0.17± 
0.03/m~ 3 of symmetric nuclear matter (SNM), at which the 
binding energy E /A per nucleon reaches its minimum. The 
empirical value E/A(po) « — 16 MeV can be extracted from 
the semi-empirical mass formula or from the extrapolation of 
binding energies of heavy nuclei. In Table [Q we display the- 
oretical calculations of saturation properties of SNM for the 
different approaches illustrated above. 

In the case of the BHF scheme with the Urbana three-body 
forces (first line of Table I), the reported values for the sat- 
uration point have been obtained by a polynomial fit to the 
numerically calculated EoS. The fit was performed with the 
constraint to reproduce the empirical saturation density of 
0.16 fm~ 3 while the corresponding energy was tuned to give 
the optimal value within a recent Energy Density Functional 
(EDF) approach based on this BHF EoS [29]. The shifts in the 
values for the density and energy at saturation from the origi- 
nal numerical values (0.17 fm~ 3 and —15.2 MeV ) have only 
marginal effects on the other physical parameters reported in 
Table I. For the BHF case with microscopic three-body forces 
(second and third lines), no application to an EDF scheme was 
performed, so the reported saturation point corresponds to the 
numerical calculations. 

In the case of the considered variational method (fourth 
line), in the quoted reference the original numerical EoS was 
corrected around saturation with an extra binding energy in 
functional form in order to reproduce the empirical saturation 
point, as reported in Table I. This correction is large as much 
as 4 MeV. We do not know how much this correction could af- 
fect the values of the other physical parameters reported in the 
table, but we used the corrected EoS, since it is the one rec- 
ommended in the original paper fl27ll and widely used in the 
literature. In particular it was used to construct the functional 
NRAPR 13011 and in many astrophysical studies on Neutron 
Stars, see e.g. ref. ifJIIl . Finally, in the DBHF approach there 
is no need to introduce three-body forces, although the satu- 
ration point is slightly shifted to higher density. We therefore 
report the calculated values without fine tuning or corrections. 

It has to be stressed that in all cases the EoS above satura- 
tion was in no way corrected. 

In the third column of the same table we show the values 
predicted for the incompressibility at the saturation point. For 
practical reasons, it is customary among nuclear physicists to 



use the definition 



d 2 E/A 



(18) 



Po 



which has the dimension of an energy. The values of the 
incompressibility extracted from the EoS in the many-body 
methods are displayed in Table Q] The values are typical of 
a pretty soft EoS at saturation, in agreement with the val- 
ues extracted from the phenomenology of monopole oscilla- 
tions ll32ll . One finds indeed that a correlation exists between 
incompressibility and position of the monopole Giant Reso- 
nance, so that, in principle it is possible to extract from the 
experimental data the value of the incompressibility in nuclear 
matter. At present, the constraints on the value of the nuclear 
matter incompressibility from the monopole excitation are not 
so tight. It can be approximately constrained between 210 and 
250 MeV (4J], though a more refined value can be expected to 
come out in the near future from additional analysis of phe- 
nomenological data. 

These values are compatible with the ones obtained from 
the experiments on sub-threshold kaon production, noticeably 
the ones from the KaoS Q and FOPI Q collaborations. 
The optimal energy for this type of investigation is close or 
even below two-body threshold, since then the only way to 
produce the kaons is by compression of the matter. Since 
at threshold the production rate increases steeply, there is a 
strong sensitivity to the value of the maximum density reached 
during the collision, and this is an ideal situation for studying 
the EoS and its incompressibility. The comparison of the sim- 
ulations with the experimental data on K + production points 
in the direction of a soft EoS. However it has to be kept in 
mind that, in the simulations, kaon production occurs at den- 
sity p > 2 — 3po, and therefore this set cannot be directly com- 
pared to the one of the monopole oscillations. In any case, 
a stiff EoS above saturation seems to be excluded from this 
analysis, as it is apparent in ref. 03511 . 

Further important properties characterizing the nuclear 
matter are the symmetry energy at saturation S(po) and its 
density derivative L, which are displayed in columns 4 and 5 
of Table J] For both those quantities experimental data do ex- 
ist, and will be discussed in the next subsections. The usual 
way of calculating them is to assume asymmetric nuclear mat- 
ter, where the proton number N p is different from the neutron 
number A^„, with N — N„ +N P . In this case, the EoS of nuclear 
matter has to be generalized. Defining 



J8 = 



N„-N p _ p„-Pp 



Nn+Np p 

as the "asymmetry" parameter, one easily gets 



|(p,/3) = j(p,0)+E sym (p)P 2 



(19) 



(20) 
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where e = (p — po) / po, Kq is the incompressibility at the satu- 
ration point, Sq = S(po) is the symmetry energy coefficient at 
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BHF.Av^ + UVIX 


0.16 


-15.98 


212.4 


31.9 


52.9 


BHF, Av 18 + micro TBF 


0.2 


-15.5 


236. 


31.3 


82.7 


BHF, Bonn B + micro TBF 


0.17 


-16. 


254. 


30.3 


59.2 


APR,Av u +UVIX 


0.16 


-16. 


247.3 


33.9 


53.8 


DBHF, Bonn A 


0.18 


-16.15 


230. 


34.4 


69.4 



TABLE I: Calculated properties of symmetric nuclear matter. 
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FIG. 1 : (Color online) Pressure as a function of baryon density for 
symmetric nuclear matter. See text for details. 
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FIG. 2: (Color online) The symmetry energy is displayed as a func- 
tion of the nucleon density. The purple zone represents the recent 
data by P. Danielewicz, whereas the different curves are the results 
of the microscopic many-body methods. 



saturation and the parameters L and K sym in the square bracket 
characterize the density dependence of the symmetry energy 
around the saturation point. The values in Table I of these 
physical parameters are obtained, for each one of the indicated 
microscopic many-body EoS, from their analytical form or by 
a polynomial fit to the numerical values reported in the origi- 
nal papers. 

B. Heavy ions phenomenology 

In the last two decades intensive studies of heavy ion re- 
actions at energies ranging from few tens to several hundreds 
MeV per nucleon (hereafter indicated as MeV/A) have been 
performed. The main goal has been the extraction from the 
data of the gross properties of the nuclear EoS. It can be ex- 
pected that in heavy ion collisions at large enough energy nu- 
clear matter is compressed and that, at the same time, the two 
partners of the collisions produce flows of matter. In princi- 
ple the dynamics of the collisions should be connected with 
the properties of the nuclear medium EoS and its viscosity. In 
the so called "multifragmentation" regime, after the collision 
numerous nucleons and fragments of different sizes are emit- 
ted, and the transverse flow, which is strongly affected by the 



matter compression during the collision, can be measured. 

Based on numerical simulations, in reference 113 611 a phe- 
nomenological range of densities was proposed where any 
reasonable EoS for symmetric nuclear matter should pass 
through in the pressure vs. density plane. The plot is repro- 
duced in Fig[T] where a comparison with the microscopic cal- 
culations discussed in Section [II] is made. The green dashed 
box represents the results of the numerical simulations of the 
experimental data discussed in ref.[36], and the brown filled 
region represents the experimental data on kaons production 
ll35ll . We notice that the EoS calculated with the BHF and the 
variational methods including UVIX three-body forces look 
in agreement with the data in the full density range. On the 
other hand, the BHF EoS obtained using microscopic TBF's 
are only marginally compatible with the experimental data, as 
well as the DBHF EoS, showing that they are too repulsive 
already at density p > 3po if the Bonn potentials are used. 
Though, it has to be stressed that all EoS are compatible with 
the data around the saturation density, i.e. their incompress- 
ibility is as soft as required by the data. However, the values of 
the incompressibility do not characterize completely the EoS, 
since it is density dependent, but in any case the analysis indi- 
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cates that the EoS al low density must be soft. 

A further constraint on the EoS is given by the symmetry 
energy, which has been extensively studied both from the the- 
oretical and experimental point of view in ref.[37J. The sym- 
metry energy is displayed in Figj2]as a function of the nucleon 
density. The purple region is the result of a recent analysis per- 
formed by P. Danielewicz on the isobaric analog states (IAS) 
in nuclei IB8I1 . This stems from the charge independence of 
nuclear interactions, i.e. strong interactions between nucleons 
in the same state do not depend on whether the nucleons are 
protons or neutrons. Therefore the energy difference between 
the ground state of a nucleus with N > Z and the isobaric 
analogs of the ground states of neighboring isobars are given 
by the symmetry energy, and the Coulomb contributions to the 
binding energy can be determined using the IAS. Many such 
states have been identified , and by fitting the available data on 
the IAS, Danielewicz and Lee obtained the constraint shown 
as a purple region in Figj2] We observe that all EoS give re- 
sults in very good agreement with the experimental data, ex- 
cept the ones of BHF with microscopic TBF at densities below 
the saturation density. 

In Figf3]we display the slope L as a function of the symme- 
try energy at saturation density, which has been widely dis- 
cussed in ref. 13711 . Several experimental data are displayed. 
The blue band represents experimental data from HIC, ob- 
tained from the neutron and proton spectra from central col- 
lisions for 124 Sn+ 124 Sn and 112 Sn+ 112 Sn reactions at 50 
MeV/A i39ll . At the same incident energy, isospin diffusion 
was investigated. We remind that isospin diffusion in HIC de- 
pends on the different N/Z asymmetry of the involved projec- 
tiles and targets, hence it is used to probe the symmetry energy 
The full red circle shows the results from isospin 
diffusion observables measured for collisions at a lower beam 
energy of 35 MeV per nucleon |j42ll . 

Transverse collective flows of hydrogen and helium iso- 
topes as well as intermediate mass fragments with Z < 9 
have also been measured at incident energy of 35 MeV/A in 
™Zn +™ Zn , 64 Zn + 64 Zn, 64 Ni + 64 Ni reactions and com- 
pared to transport calculations. The analysis yielded values 
denoted by the full squares j43ll . 

The box labelled by FRDM (finite-range droplet model) 
represents a refinement of the droplet model 044[, and includes 
microscopic "shell" effects and the extra binding associated 
with N = Z nuclei. The FRDM reproduces nuclear binding 
energies of known nuclei within 0.1%, and allows determina- 
tion of both So = 32.5 ± 0.5 MeV and L = 70 ± 15 MeV. 

In Figf3] the other boxes represent experimental data ob- 
tained from measurements of the neutron skin thickness. In 
light nuclei with N ~ Z , the neutrons and protons have sim- 
ilar density distributions. With increasing the neutron num- 
ber N, the radius of the neutron density distribution becomes 
larger than that of the protons, reflecting the pressure of the 
symmetry energy. The measurement of the neutron skin thick- 
ness is made on the stable nucleus 208 Pb, which has a closed 
neutron shell with A/ = 126 and a closed proton shell with 
Z = 82, hence it is very asymmetric and the neutron skin is 
very thick. The possibility of measurements of the neutron 
radius in 208 Pb by the experiment PREX at Jefferson Labora- 
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FIG. 3: (Color online) The derivative of the symmetry energy L is 
shown as a function of the symmetry energy at saturation Sq. See 
text for details. 

tory has been widely discussed B45I1 . The experiment should 
extract the value of the neutron radius in 208 Pb from parity- 
violating electron scattering. However, the experimental sig- 
nature is very small, and the extracted thickness has a large 
statistical uncertainty. In the next few years, a second experi- 
mental run for PREX could reduce this large uncertainty [46] . 

Recent experimental data obtained by Zenihiro et al. 04711 on 
the neutron skin thickness of 208 Pb deduced a value 8R np = 
0.21 l^o 063 F rom tne experiments constraints on the sym- 
metry energy were derived, and these are plotted in Figf3]as 
the short-dashed blue rectangular box labelled Pb(p,p). 

Last, we mention the experimental data on the Pygmy 
Dipole Resonance (PDR) in very neutron-rich nuclei such as 
68 Ni and 132 Sn, which peaks at excitation energies well be- 
low the Giant Dipole Resonance (GDR), and exhausts about 
5% of the energy-weighted sum rule 114 811 . In many models it 
has been found that this percentage is linearly dependent on 
the slope L of the symmetry energy. Carbone et al. ll49l ex- 
tracted a value of L = 64.8 ± 15.7 MeV, and So = 32.2 ±1.3 
MeV using various models which connect L with the neutron 
skin thickness. Those constraints are shown as a long-dashed 
rectangle in Fig|3]with the label PDR. 

The predictions of the different EoS are also reported in 
Figf3]as full symbols. They are distributed within a large re- 
gion and they span a wide interval in the values of the parame- 
ter L. However, the various phenomenological data are at best 
marginally compatible, and it is difficult to put well definite 
constraints on the EoS. Tentatively, from these data one can 
restrict the possible values of the symmetry energy at satura- 
tion in a limited interval, approximately 30 < So < 35 MeV, 
where all the considered EoS are actually falling. 
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FIG. 4: (Color online) The mass-radius (left panel) and the mass-central density (right panel) relations are plotted for the EoS's discussed. 
Boxes are boundaries extracted from observations, see ref. (Hf]. 




FIG. 5: (Color online) Pressure as a function of the mass-energy density in neutron star matter. The shaded areas are taken from ref. lBill . See 
text for details. 



C. Astrophysics 



A neutron star is bound by gravity, and it is kept in hydro- 
static equilibrium only by the pressure produced by the com- 
pressed nuclear matter. It is then apparent that the nuclear 
matter EoS is the main medium property that is relevant in this 
case, as can be seen in the celebrated Tolman-Oppenheimer- 



Volkoff |[30h equations, valid for spherically symmetric NS 
em /„ f\ / 47r7V 3 \ 2Gm 



dP 
dr 
dm 

~d7 



= —G- 



4nr 2 e 



(! -) 



-l 



(21) 



where G is the gravitational constant, P the pressure, £ the 
energy density, and r the (relativistic) radius coordinate. To 
close the equations we need the relation between pressure and 
density, P = P(e), i.e. just the EoS. Integrating these equa- 
tions one gets the mass and radius of the star for each central 
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density. Typical values are 1-2 solar masses (M ) and about 
10 Km, respectively. This indicates the extremely high density 
of the object. It turns out that the mass of the NS has a maxi- 
mum value as a function of radius (or central density), above 
which the star is unstable against collapse to a black hole. The 
value of the maximum mass depends on the nuclear EoS, so 
that the observation of a mass higher than the maximum one 
allowed by a given EoS simply rules out that EoS. The consid- 
ered microscopic EoS's are compatible with the largest mass 
observed up to now, that is close to 1 .97 ±0.04 M 115 111 . This 
is clearly shown in Figj4] where the mass-radius (left panel) 
and mass-central density relations (right panel) are plotted for 
all the considered EoS's as thick lines. It looks unlikely that 
this value is indeed the largest possible NS mass, and there- 
fore future observational data on NS masses could overcome 
this limit and strongly constrain the nuclear EoS. 

It would be of course desirable to have some phenomeno- 
logical data also on the radius of NS. Unfortunately this is 
quite difficult, but some tentative analysis looks promising 
152,15311 . In Fig|4]a sample of observational data taken from 
ref. ll52ll is displayed by closed thin lines for different sources, 
measured in quiescence and from thermonuclear bursts. It 
turns out that the current measurements are consistent with 
radii in the range 8-12 km and disfavor neutron stars with R 
~15 km. Those measurements are consistent with the recent 
observation of the neutron star in SAX J1748. 9-2021, which 
points to the neutron star radius in the 8-11 km range [53]. 

Additional tentative constraints on the nuclear EoS were 
obtained in a recent analysis of the data on six NS based 
on Bayesian statistical framework 115411 . Depending on the 
hypothesis made on the structure of the NS, the results are 
slightly different, as shown in Fig|5] where the quantity r p h 
is the photosphere radius. In the let panel r p /, is compara- 
ble to the neutron star radius R, whereas in the right panel 
a substantial expansion of the photosphere during an X-ray 
burst is assumed to occur. The overall allowed region where 
the EoS's should lie is displayed in Fig|5]as bounding boxes, 
where the theoretical EoS's just discussed are also reported 
as thick lines. Among the different EoS's, only the one cal- 
culated with BHF and phenomenological Urbana model ap- 
pear to be compatible with the extracted observational con- 
straints over the whole density range. It turns out that other 
microscopic EoS do not show the same agreement, in partic- 
ular the EoS with BHF and microscopic TBF and the DBHF 
EoS look too repulsive at high density. These boundaries ob- 
tained from astrophysical data are complementary to the ones 
obtained from heavy ion reactions, and illustrated in the previ- 
ous subsection. In fact, in heavy ion collisions the tested mat- 
ter is essentially symmetric, while in NS the matter is highly 
asymmetric. Considered together, the two types of constraints 
probe the density dependence of the symmetry energy. 

In relation to the high density region of the nuclear EoS, an 
additional test is on the speed of sound c$, that is required to 
be smaller than the speed of light c (causality condition). The 
speed of sound is directly connected with the incompressibil- 
ity and the energy density, according to the relativistic expres- 
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FIG. 6: (Color online) The speed of sound is plotted as function of 
the nucleon density for the EoS's discussed in the text. The dots mark 
the central density for the maximum mass of a neutron star. 
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FIG. 7: (Color online) The gravitational mass is plotted as a function 
of the baryon mass for the EoS's discussed in the text. The boxes 
indicate the boundaries coming from the simulations. 



sion 



(22) 



As such, it depends on the matter energy density and asym- 
metry. In Fig|6]we plot the speed of sound in units of c as 
a function of the density of the NS matter, according to each 
EoS. In the figure the full points on each curve indicate the 
central density of the calculated NS maximum mass for the 
given EoS. At large enough density most of the EoS's show a 
superluminal speed of sound. One finds that also the DBHF 
EoS shows a superluminal behavior. This is not completely 
surprising, since the DBHF approach is actually based on the 
three-dimensional reduction of the original four-dimensional 
Bethe-Salpeter equation, and therefore it is not fully relativis- 
tic. 

A further additional constraint on the neutron star EoS is 
provided by the observation of the double pulsar J0737-3039, 
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and the interpretation given by Podsiadlowski ll55ll . In fact, 
the gravitational mass of Pulsar B is very precisely known 
Mg = 1.249 ± 0.001 M , whereas estimates of the baryonic 
mass depend upon its detailed mode of formation. As mod- 
elled by Podsiadlowski et al., if the pulsar B was formed from 
a white dwarf with an O-Ne-Mg core in an electron-capture 
supernova, assuming no or negligible loss of baryonic mass 
during the collapse, the newly born neutron star will have 
the same baryonic mass as the precollapse core of the pro- 
genitor star, i.e. Mb — 1.366 — 1.375M . This result is dis- 
played in Fig|7]as a black box. Though, taking into account 
the uncertainty in the EoS and the small mass loss during the 
collapse, Kitaura et al. 15611 made another simulation which 
gave M B = 1.360 ± 0.002 M e , which is shown in Fig|7]by 
the dashed red box. We have calculated for each neutron star 
matter EoS the relation between the gravitational and bary- 
onic mass, and these are displayed in Fig|7]by the straight 
curves. We notice that the results of all microscopic EoS's 
agree very well with the result of Podsiadlowski, at variance 
with the calculations based on the phenomenological Skyrme 
forces discussed in ref.(|4|, where agreement was found with 
the result of Kitaura et al. f56Tl . which assumed small mass 
loss during the collapse. 

IV. CONCLUSIONS 

We have presented a systematic confrontation of the nu- 
clear Equation of State, obtained within different microscopic 
many-body methods, with the available constraints coming 
from phenomenology. The latter are extracted from labora- 
tory experiments as well as from astrophysical observations. 
Both nuclear structure and heavy ion collisions data were con- 
sidered, along the same lines of the analysis on the Skyrme 
forces reported in ref. [4]. Astrophysical observational data 
included the measures of NS masses, some hints on the radius- 
mass relation from ref. 1B2I EH and the constraints on the 
EoS presented in ref. ll54ll . obtained from the analysis of tran- 
sient phenomena in six NS's. Some theoretical constraints, as 
the requirement of a sub-luminal speed of sound, were also 



considered. If one takes literally all the constraints, among 
the considered microscopic EoS only one passes all the tests. 
However, these phenomenological constraints are affected by 
uncertainties, which are difficult to estimate quantitatively on 
a firm basis. The conclusion that one can draw from this anal- 
ysis is twofold. Firstly, despite the differences between the 
considered microscopic EoS, their overall predictions do not 
show major discrepancies with one another as well as with 
the phenomenological constraints. In other words, the calcu- 
lated microscopic nuclear EoS, based on different many-body 
methods, are not in contradiction with the phenomenological 
constraints. This is not at all an obvious result. Secondly, 
from the analysis it appears that one can explain reasonably 
well all the data with a microscopic EoS that includes only 
nucleonic degrees of freedom, in particular no exotic compo- 
nents in NS are needed. On one hand, it is likely that this 
conclusion will be disproved in the future with the expected 
new data, in particular on the NS maximum mass. Indeed, if 
the evidence for a NS mass of 2.5 solar mass [57] will be con- 
firmed, this will rule out most of the considered microscopic 
EoS and will introduce a serious and fundamental issue in the 
physics of NS. On the other hand, it has been shown that ex- 
otic matter like hyperons or quarks should appear in NS 1I22I . 
which will strongly affect the EoS and challenge the theory of 
high density nuclear matter. In the future one can expect that 
the interplay between theory and observations will continue to 
play a major role in the worldwide effort of determining the 
nuclear EoS. 
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